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Abstract 

We reconsider the scenario in which the knee in the cosmic ray spectrum 
is explained as due to a change in the escape mechanism of cosmic rays from 
the Galaxy from one dominated by transverse diffusion to one dominated by 
drifts. We solve the diffusion equations adopting realistic galactic field mod- 
els and using diffusion coefficients appropriate for strong turbulence (with a 
Kolmogorov spectrum of fluctuations) and consistent with the assumed mag- 
netic fields. We show that properly taking into account these effects leads to 
a natural explanation of the knee in the spectrum, and a transition towards 
a heavier composition above the knee is predicted. 



1 Introduction 

A well established feature of the full cosmic ray (CR) energy spectrum is that it has 
a power-law behavior with a steepening taking place at the so-called knee, corre- 
sponding to an energy E^ nee ~ 3 x 10 15 eV. Although the knee of the CR spectrum 
has been known since more than four decades, none of the numerous models pro- 
posed so far to explain this feature has managed to become broadly accepted. Some 
proposals focus on a possible crossover between different acceleration mechanisms 
below and above the knee [1, 2, 3], or exploit the possibility of a change in the 
particle acceleration efficiency [4, 5, 6]. Other hypotheses include the nuclear pho- 
todisintegration at the sources [7, 8], the recent explosion of a single source [9], and 
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leakage from the Galaxy due to a change in the confinement efficiency of CRs by 
galactic magnetic fields [10, 11, 12]. Among the latter, Ptuskin et al. [12] in partic- 
ular consider the knee as due to a crossover from a diffusive regime dominated by 
transverse diffusion at low energies to another dominated by drifts, i.e. by the Hall 
(antisymmetric) diffusion, above the knee. Since the diffusion coefficients determine 
the residence time of the CRs inside the Galaxy, this means that the observed slope 
of the CR energy spectrum will differ from the slope of the original average source 
spectrum just due to the energy dependence of the diffusion coefficients. Hence, the 
scenario just mentioned naturally accounts for a change in the spectral slope, arising 
from the diverse energy dependence of the two different diffusion coefficients. 

In Ref. [12] some important simplifying assumptions had to be done to make the 
differential equations determining the CR densities more tractable. In particular, 
a simplistic spatial dependence of the regular field and of the diffusion coefficients 
was adopted, but it was however found that the spectral slope beyond the knee 
was quite sensitive to the magnetic field configuration adopted and to the spatial 
distribution of the sources. Also, the expressions for the diffusion coefficients used 
were actually only valid in the limit of small turbulence. Anyhow, a feature similar 
to the knee was obtained, but to reproduce the correct value of Ek nee a magnetic field 
spectrum flatter than Kolmogorov was needed. This last requirement resulted from 
the fact that the perpendicular diffusion coefficient (the only one relevant for the 
determination of the CR densities below the knee) was normalized at low energies 
(few GeV) to the average value resulting from a comparison with the results of a 
simplified leaky box model of the Galaxy, and was then extrapolated up to E^nee 
using the fact that its energy dependence is directly related to the spectral index of 
the random magnetic fields. 

In the present work we want to further elaborate on the proposal of Ref. [12], 
but using more realistic configurations for the galactic magnetic fields (including e.g. 
field reversals, which enhance field gradients and can hence affect the CR drifts). 
We also consider diffusion coefficients appropriate for strongly turbulent magnetic 
fields, since strong turbulence from disturbances in the interstellar plasma is thought 
to be the actual picture of magnetic fields in the Galaxy, and we will directly adopt 
diffusion coefficients normalized from the results of numerical simulations valid in the 
relevant energy range of the knee and consistent with the magnetic fields adopted. 

The beauty of this scenario is that in it the presence of a break in the spectrum 
does not require any special assumption, but just follows from properly taking into 
account well established properties of the propagation of charged particles in reg- 
ular and turbulent magnetic fields. Moreover, this break naturally happens at the 
energies at which the knee is actually observed. 
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2 The diffusion of cosmic rays in the Galaxy 



2.1 Regular and random galactic magnetic fields 

The Milky Way, as well as other spiral galaxies, is known to have a large-scale reg- 
ular magnetic field B (x), which is determined by measurements of the rotation 
due to the Faraday effect of the polarization plane of radiation from pulsars and 
extra-galactic radio sources [13, 14]. The magnetic lines are observed to follow the 
spiral arms, but there is a variety of models to describe the regular magnetic field 
structures. Models of magnetic fields of spiral galaxies are distinguished between 
bisymmetric (BSS) and axisymmetric (ASS), defined so that BSS (ASS) fields are 
antisymmetric (symmetric) with respect to a n rotation. For a spiral galaxy with 
two arms, for instance, in BSS models the field reverses its sign between the arms, 
while in ASS models it keeps the same orientation. In the case of the Milky Way, the 
regular magnetic field in the galactic disk is predominantly toroidal, being aligned 
with the spiral structure and having reversals in direction between nearby arms. 
Hence, the regular field in the disk is usually described with a BSS model. Further- 
more, magnetic field models can be classified as symmetric (S) or antisymmetric 
(A) with respect to the galactic plane, and the regular disk component is believed 
to be symmetric with respect to the galactic plane. The presence of substantial 
non-thermal radio emission at high galactic latitudes points to the existence of an 
extended corona with non-negligible magnetic fields. This magnetic halo extends 
possibly up to 10 kpc from the galactic plane. Although the orientation of the reg- 
ular field in the galactic halo is not known, the toroidal component appears to be 
dominant. The magnetic field structure of the halo may be similar to that of the 
galactic disk (for instance if its origin is related to a galactic wind which somehow 
extends the field disk properties into the halo), but it may also be different since it 
could be independently generated through a dynamo mechanism. 

Taking into account the main features of the regular magnetic field in the Galaxy 
mentioned above, we will assume that the propagation region of the galactic CRs 
is a cylinder of radius R = 20 kpc and height 2H, with a halo size H = 10 kpc. 
Moreover, since the pitch angle of the spiral arms is quite small (< 10°), it is a good 
approximation to consider that the regular magnetic field is directed in the azimuthal 
direction. We will hence adopt cylindrical coordinates (r, 0, z) throughout, taking 
B = £>o0 and assuming for simplicity that the system has azimuthal symmetry. In 
passing, notice that a 0— independent azimuthal field is necessarily divergenceless. 
In order to reproduce the field reversals of the BSS-S model, we will adopt for the 
regular disk component the expression 
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where r c = 4 kpc is a core radius which smoothes out near the galactic center the 
overall 1/r behavior, — 1 kpc is the adopted magnetic disk scale height, and the 
value of B d is chosen such that BQ lsk (r obs , z obs ) = —2 fiG at our observation point 
i r obs = 8.5 kpc, z bs = 0) 1 . Notice the minus sign arising from the fact that the 
local galactic magnetic field is nearly directed towards the longitude / ~ 90°, i.e. in 
the —(f) direction. This regular field component has reversals for r = 4, 8, 12 and 
16 kpc. For the regular halo component, we will consider it to be described either 
by a field structure with radial reversals, namely 

B^ lo (r,z)=B h , 1 -= 1 sin ( ^— W ( ^— ) , (2) 

yfl + (r/r c ) 2 ft + (z/z h ) 2 V 4k P c / V lk P c / 

or by a model without radial reversals, as given by 

where Zh = 5 kpc is the adopted magnetic halo scale height, and Bh is such that 
BQ alo (r obs , z obs ) = ±1 A*G. The function 1Z is just unity for the symmetric halo case, 
and it is taken as TZ(z) = sign(^) in the case of a halo model antisymmetric with 
respect to the galactic plane. 

Eqs. (l)-(3) provide simple but yet realistic galactic magnetic field models, con- 
taining the characteristic features (such as the reversals) which we believe are im- 
portant for the determination of the drift effects. 

Superimposed to the regular magnetic field structures of the galactic disk and the 
extended halo, an irregular magnetic field B r (x) due to turbulence in the interstellar 
plasma is known to exist, with the largest eddies having a scale L max ~ 100 pc. 
This random field component has a strength comparable to the field of the regular 
component, namely (|B r |)/ (|B |) ~ 1-3 locally. The existing observational data [15] 
show that the spectrum of inhomogeneities may be the same for the density of the gas 
and for the magnetic field, and that it is close to a Kolmogorov spectrum. Hence, 
the random magnetic field can be assumed to be given by a power-law spectrum 
with Fourier components giving rise to a magnetic energy density dE/dk oc k~ 5 / 3 , 
for k > 27r/L max . The random field intensity will be considered to possess a smooth 
spatial dependence analogous to the envelope of the regular fields assumed above, 
namely 

\B r (r,z)\=B* r - 1 - T ^= , (4) 



1 Thc th 2 (r) factor multiplying the regular magnetic field components just insures that there 
are no singularities in the r — axis, something which would otherwise artificially introduce a 
singular drift current along that axis (see below). 
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where z r = 3-5 kpc and B* is such that |B r (r 6 S , z b s )\ = 2-4 fiG. Since we are 
concerned with highly relativistic particles of speed v ~ c, which is much larger 
than the Alfven velocity va characteristic of the magnetic wave propagation, it is 
well justified to assume the magnetic disturbances to be static, and thus to neglect 
electric fields. 

2.2 The diffusion equation 

On scales larger than the scale of magnetic irregularities, the cosmic ray transport 
across the turbulent galactic magnetic fields can be well described by means of a 
diffusion equation, in which the components of the diffusion tensor depend on both 
the regular field and the random turbulence. Neglecting nuclear fragmentation and 
any energy losses, the steady-state diffusion equation for the cosmic ray density 
iV(x) reads 

V-J = Q(x) , (5) 
where Q is a source term and the CR current is given by 

J, = -A,(x)V,iV(x) . (6) 

In general, the diffusion tensor Dy can be written as 

Dij = (l>\\ - £>_)_) bibj + D ± 5ij + D A e ijk bk , (7) 

where b = B /|B | is a unit vector along the regular magnetic field, 5ij is the Kro- 
necker delta symbol, and is the Levi-Civita fully antisymmetric tensor. The 
symmetric terms contain the diffusion coefficients parallel and perpendicular to the 
mean field, D\\ and D±, which describe diffusion due to small-scale turbulent fluc- 
tuations. The antisymmetric term contains the Hall diffusion coefficient Da, which 
is responsible for the macroscopic drift currents. The diffusion equation can be 
rewritten as 

-ViDsijVjN + UiViN = Q , (8) 

where Dsij = — D±^J bibj + D±5ij contains the symmetric part of the diffusion 
tensor, and U; L = — e^Vfe (-Da^j) is called the drift velocity. Notice also that the CR 
current in Eq. (6) may be written as 

J = -D ± V±N -D\\V\\N + D A b x ViV , (9) 

where Vy = b(b • V) and Vj_ = V — Vy. Hence, the macroscopic drift current 
is orthogonal to the overall direction of the regular field and to the CR density 
gradient. 
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Under the assumption mentioned previously that the regular magnetic field is 
directed in the azimuthal direction (i.e. b r = b z = 0, with = ±1), and that the 
system has azimuthal symmetry {d/d<p = 0), Eq. (5) simplifies to 




) 



l( D ^ +u 'l + u 4] N{r > z) = Q{r ' z) • (10) 



where the drift velocities are 



u r = 



d{D A h) 

dz 



(11) 



and 



1 d{rD A b^) 



(12) 



r dr 



These expressions for the drift velocities already make apparent the possible rele- 
vance of including field reversals in the regular fields, since this certainly affects the 
spatial gradients of the diffusion coefficients. 

At the surface E corresponding to the boundaries of the halo, the CR density has 
to become negligible, since particles can escape freely into the intergalactic space. 
The boundary conditions that apply to Eq. (10) are thus N\^ = 0. 

Notice that the diffusion treatment of the CR propagation will be adequate as 
long as the Larmor radius, which for a relativistic particle of energy E and charge 
Z is given by 



is smaller than the typical scale height of the magnetic field. Hence, for CRs trapped 
in the disk of the Galaxy, this requires E/Z < few x 10 17 eV, and somewhat larger 
values could result if the magnetic fields in the halo were sizeable. Clearly the knee 
region is well inside this domain. 

2.3 Calculation of transverse and Hall diffusion coefficients 

The components of the diffusion tensor that define the CR diffusion process through 
Eqs. (10)-(12) depend on both the regular field and the random turbulence, but un- 
fortunately a comprehensive theory capable of determining the diffusion coefficients 
from our knowledge about the magnetic fields is lacking. 

From quite general considerations concerning the ensemble average of the two- 
time particle velocity correlations, i.e. Rij(t) = (vi(t )vj(t + t)), applied to the 
Taylor-Green-Kubo formula [16], the following theoretical relations involving the 
diffusion coefficients are expected to hold: 
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vr L ut ± 

D± = -YT(^f (15) 

and 

vtl {ur A f 
3 l + (ujT A f ^ ] 

where v ~ c is the particle speed and uj = i>/t\l is the angular gyrofrequency 2 . 
It should be noticed that three different timescales tn , r± and appear here as 
the effective timescales for the decorrelation of particle trajectories. Indeed, the 
trajectory decorrelation due to a given component of the diffusion tensor will be in 
general associated to some particular physical processes which may not always be 
contributing to the other components of the diffusion tensor. For instance, for low 
turbulence levels the field line random walk appears as the dominant contribution 
to the transverse diffusion, but it plays no major role in the parallel diffusion. 

There is however no general theory providing the decorrelation timescales. In 
particular, this turns out to be a severe restriction for the study of particle diffusion 
under highly turbulent conditions. It should be noticed that in treating CR diffu- 
sion in the Galaxy we are precisely concerned with high turbulence levels, not only 
because the mean random magnetic field is of the order of the mean regular field, 
but also because when there are reversals in the orientation of the regular field, there 
exist regions in which the regular field becomes negligible and hence inside them the 
turbulence largely prevails. 

To go beyond the general expressions in Eqs. (14)-(16) requires then to make 
some simplifying assumptions about the timescales for decorrelations. For instance, 
in Ref. [16] a common decorrelation timescale r was assumed to hold for both the 
perpendicular and the antisymmetric diffusion (r = r± = ta)- It was then assumed 
that the major contribution to perpendicular diffusion is due to field line random 
walk and, based on simple considerations, the Ansatz 

2 7V 

^=o7T- (17) 

was adopted to relate the decorrelation timescale to the diffusion coefficient Df trw 
associated with the field line random walk. Provided either an analytical or numer- 
ical estimate for Dfi rw , which depends on both the form of the turbulence spectrum 
and the turbulence level assumed, D± and D A are thus obtained. However, the 
numerical study of Ref. [17] shows that the Ansatz proposed in [16] leads to results 
clearly inconsistent with the simulations, so it may be concluded that Eq. (17) is 
not appropriate to model all effects contributing to the randomization of the particle 
trajectories. 

2 Thc sign of Da in Eq. (16) corresponds to the diffusion of positively charged particles, and 
the opposite sign will hold for negatively charged ones. 
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Expressions similar to Eqs. (14)-(16), but in which only a single timescale r 
appears for the three diffusion coefficients, were obtained in several previous ana- 
lytic studies that assumed a single scattering process to be responsible for all the 
decorrelations. This happens, for instance, by considering a nearly isotropic particle 
distribution under an homogeneous scattering process, in which r represents the 
characteristic time of isotropization [18]. A similar structure was also found for col- 
lisional diffusion coefficients in a plasma in thermal equilibrium [19]. Furthermore, 
regarding r as the scattering time, these expressions can also be obtained by treating 
the magnetic disturbances as hard-sphere scattering centers [20], and is known as 
the "classical scattering result" [21, 22, 23]. 

Very recently, Casse et al. [24] presented extensive Monte Carlo measurements 
of parallel and transverse diffusion of relativistic test particles propagating in strong 
stochastic Kolmogorov fields, extending previous numerical results [17]. From these 
investigations, it turns clearly out that the extrapolation of the previously mentioned 
theoretical models to account for particle diffusion in strongly turbulent fields is 
generally not applicable. For instance, the numerical simulations [17, 24] show that, 
while Eq. (14) for parallel diffusion appears to agree well with the obtained results, 
Eq. (15) for transverse diffusion (with r = t\\ = r±) clearly fails, as expected from 
the fact that field line random walk is in general not negligible. 

Given these facts, what we will do instead is to directly extract the transverse 
diffusion coefficient from the numerical results of Ref. [24]. These results will be 
applicable in the regime of interest for the study of the knee region, since the relevant 
energy parameter, p = 2nri/ L max , takes the value p ~ 0.1 for protons near the knee, 
which is just inside the range [0.001,10] explored in Ref. [24]. 

In Ref. [24] D± was obtained as a function of rigidity and turbulence level (their 
Figure 5), were the turbulence level is defined as r\ = B^/(B^ + Bq). It turns out 
that an approximate expression for D±, valid in the range 0.1 < i] < 1, is just 



where ri is an effective Larmor radius defined as in Eq. (13) but replacing |B 



y Bq + B%. It should be remarked that, although this expression is obtained from 
numerical simulations, the dependence on rigidity through a power 1/3 appears 
also in theoretical expressions (valid in the limit of small turbulence) under the 
assumption of a Kolmogorov spectrum for the random fields [12]. 

Regarding the drift effects, there are at present no numerical studies concerning 
the Hall diffusion under turbulent conditions, so that some additional assumptions 
had to be made. Notice first that from Eq. (16) we can obtain the value of D A 
in the two limiting cases of very strong or weak turbulence. In the first case, i.e. 
for ujta — > 0, one has as expected that Da — > 0, so that drift effects become small 




(18) 
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(this limit is of particular interest in the regions with field reversals, in which the 
turbulence will tend to suppress the drift effects). On the other hand, in the limit 
of very weak turbulence (ur A ^> 1), the obtained expression tends to the well- 
known result Da = vtl/3 of the quasi-linear approximation. In order to obtain 
a smooth transition between the two limiting regimes, and lacking theoretical or 
numerical results for ta, we will just adopt Eq. (16) for Da but assuming here 
a common timescale for all decorrelation processes. Under the assumption of a 
common decorrelation timescale, it turns out in particular that 

^ = l + (^) 2 , (19) 

and this ratio was obtained as a function of rigidity and turbulence level in Ref. [24] 
(see their Figure 6). The ratio D\\/D± turns out to be almost independent of the 
energy in the range studied, and using Eq. (19) we find that the simple approxi- 
mate relation ljt ~ 4.5(|B |/|B r |) L5 reasonably reproduces the numerical results 
for Dn/D±. This result for lot may now be replaced in Eq. (16) (after changing 
ta — > t) and this provides an expression for the Hall diffusion coefficient. The re- 
sulting expression for D A has the proper behavior in the limits of very strong or very 
weak turbulence, and the criterium adopted to obtain cut just allows to get a smooth 
transition between the two regimes. We do not expect the general results obtained 
to change qualitatively if a different prescription to regulate D A were adopted. Any- 
how, further investigations concerning the Hall diffusion under turbulent conditions 
are currently under progress. 



3 Results and comparison with observations 

We have assumed the regular galactic magnetic field to be composed of a disk 
component with radial reversals and symmetric with respect to the galactic plane 
(given by Eq. (1)), and a halo component which could be described either by a 
field structure with radial reversals (Eq. (2)) or by a model without radial reversals 
(Eq. (3)). In the former case, it is supposed that the halo is an extension of the disk 
magnetic field, and hence it will be considered to be also symmetric with respect 
to the galactic plane and directed parallel to the disk component throughout the 
Galaxy. In contrast, in the latter case all relevant variants concerning symmetry and 
orientation will be explored, since this case may correspond to an independently 
generated halo field. Hence, in this case four different configurations arise from 
choosing the orientation (through the sign of Bh) and the symmetry with respect to 
the galactic plane (through the function 1Z) . 

Adopting a certain galactic field model and assuming a given source distribution 
Q(r,z), the CR diffusion is determined by solving the diffusion equation given by 
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Figure 1: Spatial distributions of the CR density within the Galaxy, assuming a 
halo without radial reversals and symmetric with respect to the galactic plane, for 
a source localized at r s = 6 kpc. Notice that the sharp peaks occurring around the 
source location have been truncated, (a) At low energies below the knee the CR 
particle diffusion is mainly due to small-scale turbulent fluctuations (governed by 
D±). (b) At high energies above the knee, the macroscopic drift currents lead to 
significant large-scale asymmetries in the CR density distribution. 

Eqs. (10)-(12). Following Ref. [12], the sources are assumed to lie on a thin disk of 
height 2h s = 400 pc, so that the source term takes the form Q(r, z) = q(r)9(h s — \z\), 
where is the Heaviside function. Concerning the radial source distributions, both 
localized 3 (i.e., q(r) oc 5(r — r s ), with r s = 3 and 6 kpc) and extended (q(r) = const.) 
distributions have been studied. The resulting differential equation has been solved 
numerically using the Crank-Nicolson discretization and the Liebmann method [25] . 

In order to illustrate the role of drifts in the CR diffusion, let us first consider 
a localized source at r s = 6 kpc, and adopting for the random field an amplitude 
|B r (r fe s , z b s )\ = 2 fiG and a scale height z r = 3 kpc. Figures 1(a) and (b) show the 
spatial distribution of the CR density within the Galaxy, assuming a halo without 
radial reversals and symmetric with respect to the galactic plane, for energies both 
below and above the knee. As expected from the preceding discussions, it is observed 
that at low energies the CR particle diffusion is mainly due to small-scale turbulent 
fluctuations (governed by D±), while at higher energies the bulk of the CRs is 
driven according to the macroscopic drift currents, leading to significant large-scale 
asymmetries in the CR density distribution. 

The influence of the drift effects on the CR spectrum at our observation point 
depends on the galactic field model, on the particular set of parameters considered 
within the given model (i.e. the field amplitudes and scale heights assumed for both 

3 These sources actually correspond to rings around the galactic center, since a non-central 
point-like source would break the azimuthal symmetry assumed. 
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Figure 2: Plots of the normalized differential energy spectrum versus E/Z for the 
different regular galactic field models investigated. R—S denotes the halo with radial 
reversals and symmetric with respect to the galactic plane, while NR — S(A) refers 
to the halo without radial reversals and symmetric (antisymmetric) with respect 
to this plane. Furthermore, the orientation of the halo relative to the regular disk 
component is distinguished by means of the sign of Bd ■ Bh, as indicated. 

the regular and random field components), and on the source distribution. Notice 
that the relevant parameter in the diffusion coefficient is the Larmor radius of the 
CRs, implying that the impact of drift effects will depend on the CR rigidity, and 
hence on the value of E/Z. 

Figure 2 shows the normalized CR flux as a function of E/Z for the different 
regular galactic field models described in Sect. 2.1 with the same source and local 
random field amplitude as in Figure 1. The random field scale height considered is 
z r = 5 kpc. Notice that in this and the following figures, the CR flux appears mul- 
tiplied by E lz+1 / 3 , where jz is the spectral index corresponding to the component 
of charge Z at the source (see Eq. (20) below). Since transverse diffusion depends 
on rigidity through a power 1/3 (see Eq. (18)), it turns out that (f)E lz+l ^ is flat at 
low energies well below the knee. R — S denotes the halo with radial reversals and 
symmetric with respect to the galactic plane, while NR — S(A) refers to the halo 
without radial reversals and symmetric (antisymmetric) with respect to this plane. 
Furthermore, the orientation of the halo relative to the regular disk component is 
distinguished by means of the sign of B d - B h . 

It can be seen that most of the models considered show the appropriate tendency 
to account for the knee, since a steepening in the spectrum appears at the right value 
of E/Z (i.e. for E ~ E^ nee for protons, and higher for heavier nuclei). Drift effects 
are particularly enhanced for the antisymmetric halo could be anticipated 
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E/Z [eV] 

Figure 3: Plots of the normalized differential energy spectrum versus E/Z for the R— 
S case and different source distributions, namely homogeneously extended (q(r) = 
const.) and localized sources (q(r) oc S(r — r s ), with r s = 3 and 6 kpc). 

from Eq. (If) since for ~ ±sign(z) one has u r ~ ^2.0^5(2:) on the galactic 
plane, and this additional contribution to the drifts significantly affects the results. 
Moreover, the halo contribution to the vertical drift velocity u z will be converging 
(or diverging) from both sides of the galactic plane (rather than being similarly 
oriented as in the symmetric models). One obtains in this case a more pronounced 
knee for Bd ■ Bh > 0, while a flattening of the spectrum just before the knee is 
observed for Bd ■ B^ < 0. The pronounced bump observed in Figure 2 for this 
case is in disagreement with the observed spectra, so that this model would be 
observationally disfavored. However, the possible existence of a smaller bump just 
below the knee in the actual data was noticed by Kempa et al. [26] many years 
ago, although it still remains somewhat controversial. It is interesting however that 
features of this kind may arise from drift effects in some particular configurations of 
regular galactic magnetic fields. 

Figure 3 shows the R — S case for an extended source distribution q(r) = const., 
as well as for sources localized at r s = 3 kpc and r s = 6 kpc, with the same random 
field as in Figure 1. As may be expected, the drift effects appear stronger the farther 
the source is located, since drift effects can remove high energy CRs from the galactic 
plane all the way from the source to the observer. The same qualitative behavior 
has also been observed in the other cases studied. 

Figure 4 displays the effects of varying the parameters relevant to the random 
field. The data correspond to the R — S case and the random field parameters 
take the values |B r (r Q 5 S , z b s ) | = 2 or 4 fiG and z r = 3 or 5 kpc. As expected, the 
drifts are clearly stronger when the random field is suppressed, either by decreasing 
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Figure 4: Plots of the normalized differential energy spectrum versus E/Z for the 
R — S case and a source localized at r s = 6 kpc. The random field parameters take 
the values |B r (r obs , z ohs ) \ = 2 or 4 /iG and z r = 3 or 5 kpc, as indicated. 

its overall amplitude or its scale length. Again, the same qualitative behavior has 
been observed in the other relevant field configurations, and is due to the increased 
suppression of Da under stronger turbulence conditions. 

To compare our results with experimental data, it is assumed that the differential 
fluxes of different nuclei emitted by the sources have power-law distributions, i.e. 

<tP z {E) = , (20) 

where 1 < Z < 26. The intensities $^ and spectral indices 7^ were taken from the 
detailed knowledge about CR mass composition below the knee (and for energies 
per particle above ~ few Z x 10 10 eV) [27], and they were then extrapolated to 
energies above the knee. Making the convolution of the spectrum of each component 
of charge Z with the modulation effects produced by the drifts discussed previously 
and summing over all components from hydrogen to iron, we arrive at the all-particle 
CR flux given by the model, which can then be compared with the observed spectra 
measured by different experiments (KASCADE, CASABLANCA, DICE, PROTON 
satellite, and Akeno) [28, 29, 30, 31, 32]. 

Figure 5 shows the main contributions to the total CR flux coming from protons 
and nuclei of helium, oxygen and iron, calculated in the R — S case for a source lo- 
calized at r s = 6 kpc, with the random field parameters as in Figure 1. Figure 6(a) 
shows the total CR flux corresponding to this case compared with some observa- 
tions, and clearly a very satisfactory agreement is found between the predictions 
and the actual experimental data. As discussed above and shown in preceding fig- 
ures, varying model parameters and considering other galactic field models enables 
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Figure 5: Plots showing the main contributions to the total CR flux coming from 
protons and nuclei of helium, oxygen and iron, calculated in the R — S case for a 
source localized at r s = 6 kpc. 



to obtain somewhat different spectra. We do not attempt here to make a detailed 
quantitative test of the different models, both because the different data still have 
sizeable statistic (and systematic) uncertainties, and because the predictions also 
rely on several simplifying assumptions. For instance, the transport equation was 
regarded as purely diffusive, neglecting spallation of nuclei in the interstellar medium 
or convection due to the possible existence of a galactic wind. Also, the spatial de- 
pendence of the various components of the magnetic field in the Galaxy were taken 
as smooth as possible in order to simplify the integration of the differential equa- 
tions determining the CR densities. Although the adopted models are plausible and 
capture the main expected features of the galactic magnetic fields, certainly the real 
ones will be less idealized. Moreover, a (sensible) prescription for dealing with the 
Hall diffusion under turbulent conditions had to be adopted, due to the lack of ap- 
propriate numerical studies concerning this issue at present. Hence, an accurate fit 
to the experimental data is not to be expected, although it is reassuring to observe 
that the predictions within this scenario tend to exhibit a remarkable agreement 
with observations. 

Figure 6(b) shows a plot of the mean mass composition (In A) versus E corre- 
sponding to the same case as before. For comparison, we also show the CR mass 
composition measured by different experiments [28, 29, 30, 31, 32]. Since leakage 
from the Galaxy depends on E/Z and is more effective for lighter particles, this sce- 
nario can easily account for a change in the composition towards heavier nuclei above 
the knee, which is apparent on some of the existing data and in particular in the 
latest KASCADE observations [32]. However, since above 10 14 eV/nucleus the infor- 
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Figure 6: Comparison between experimental observations and model predictions 
obtained in the R — S case for a source localized at r s = 6 kpc: (a) plot of the total 
CR flux versus the energy per particle; (b) plot of the mean mass composition (InA) 
versus the energy per particle. 

mation about the CR mass composition has to be obtained indirectly from extensive 
air shower observations at ground level, still significant discrepancies are found be- 
tween the results of experiments using different techniques. Indeed, the CR mass 
composition beyond the knee is not definitively established yet, with some observa- 
tions [28] suggesting that it could turn lighter while others finding that the heavier 
components become dominant [31, 32]. In Figure 6(b) it can be seen that some data 
sets suggest that the mean mass composition has a dip starting at E ~ 10 15 eV, 
while it grows steadily towards heavier components beyond E ~ 3 x 10 15 . A dip cor- 
responding to the composition becoming lighter could possibly result from a bump 
in the spectrum of the proton component in the same energy range. The models 
predict on the other hand that for increasing energies only the heavier elements 
persist, since the lighter ones are efficiently drifted away for E/Z 3> E knee . 

As a summary, we have considered in detail the diffusion and drift effects of 
CRs in realistic configurations of the galactic magnetic fields, adopting diffusion 
coefficients appropriate for the high levels of turbulence present in the Galaxy and for 
a Kolmogorov spectrum of fluctuations. On the one hand, the larger field gradients 
associated to the field reversals present between spiral arms enhance the drift effects, 
but on the other hand a suppression of associated to the high turbulence present 
(particularly near field reversals) tends to reduce the drifts, so that the two effects 
compete among each other. We used diffusion coefficients obtained directly through 
numerical simulations valid in the relevant energies, and found that the correct 
value of E knee can easily be obtained. Although the slope of the spectrum of each 
nuclear component changes by a factor Ao ~ —2/3 from below to above ZEk nee , the 
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envelope of the total spectrum shows a milder change in slope, which is consistent 
with the observed spectrum. A progressive transition towards a heavier composition 
results in this scenario. 

Certainly a better knowledge of the CR spectrum and especially of the CR mass 
composition around and above the knee is of major importance to finally decide 
which is the correct mechanism responsible for this feature. 
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